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ABSTRACT 

Aims. The observations of TeV blazars published recently show an unexpected quadratic or even cubic correlation between the X-ray 
and gamma-ray emission. A standard model of the synchrotron self-Compton emission of a compact source inside a jet is not able 
to explain such a correlation. Therefore, we propose an alternative scenario where the emission of at least two independent compact 
components is observed at the same time. 

Methods. We compare two different models. The first model assumes the injection of relativistic particles into a downstream region 
of a shock wave inside a jet that creates the emitting source. The model precisely describes the evolution of the particle energy 
spectrum inside the source and takes into account a light-crossing time effect for the produced radiation. The second model assumes 
an intrinsically constant emission of a homogeneous source that travels inside the jet along a curved trajectory, where the activity is 
produced simply by different values of the source's Doppler factor. To verify the two models we use recentlu published observations 
ofMrk421. 

Results. Our simulations show that simultaneous radiation of at least two independent sources, where the first source dominates the 
emission in the X-ray range and the second source radiates strongly in the gamma-ray range, can explain the observed correlations. 
However, the injection model provides inadequate results because it gives different values for the correlation of the rise and decay of 
a flare. This problem is negligible in the scenario that uses the Doppler boosting effect. Therefore, this approach yields much better 
results. 

Key words. Radiation mechanisms: non-thermal - Galaxies: active - BL Lacertae objects: individual: Mrk 421 



1. Introduction 

Activity of TeV blazars is usually observed simultaneously in 
the X-ray and the gamma-ray range (e.g. Catanese et al. 1 1 9971 
Pian et al. [19981 S ambru na et al. |2000l Tak ahashi et al. 120001 
Krawczynski et al. 2002 Donnarumma et al. 120091 , Therefore, 
the activity can be analysed by a simple comparison of the two 
different light curves. Such a comparison has shown in a few 
cases a surprisingly precise correlation between the evolution of 
the X-ray and the gamma-ray emission. This provides an excel- 
lent opportunity to test theoretical models for very high energy 
emission of blazars. 

To describe the correlation is convenient to assume that the 
evolution of the X-ray or gamma-ray flux during an outburst 
can be approximated by a power-law function (Fx °c t s and 
F Tev K t c ) with different values of s, c for the rise and decay 

1 / s 

phase of the activity. Rewriting the first relationship as t oc F x 
and substituting this to the second proportionality we can define 
the correlation 



oc F 



cj s 



(1) 



which is also a power-law function with the index x = c/s. Such 
a simple definition can describe precisely the correlations ob- 
served in TeV blazars. 

The analysis of the light curves produced by different blazars 
shows that there is no unique value of the correlation slope. The 
activity of Mrk 501 observed in April 1997 (Catanese et al. l 19971 1 
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gives x = 1.7 1+0. 50 for the comparison between the 50-1 50 keV 
light curve obtained by the OSSE experiment and the gamma-ray 
observations obtained above 350 GeV by the Whipple telescope. 
On the other hand the comparison between the 2-10 keV obser- 
vations made by the RXTE-ASM experiment and the gamma 
rays detected by the Whipple for the same activity of Mrk 501 
gives x = 2.69+0.56 (Katarzynski et. al. 2005). Moreover, obser- 
vations of Mrk 501 conducted by the RXTE-PCU experiment (2- 
20 keV) and the gamma-ray telescopes (HEGRA and Whipple) 
show a linear correlation (x = 0.99 + 0.01) for the data obtained 
in May 1997 and an almost quadratic relation (x = 2.07 + 0.27) 
for the observations made in June 1998 (Gliozzi et al. I2006t . 

The observations of Mrk 421 made in March 2001 (Fossati et 
al. 12008) show several outbursts observed simultaneously in the 
X-rays (RXTE) and the gamma rays (HEGRA, Whipple). The 
detailed data analysis performed by the authors of the observa- 
tions shows at least in one case more than a quadratic correlation 
(x > 2) and in two cases significantly more than linear relation 
(x > 1). 

Recent observations of PKS 2155-304 made by Chandra and 
the H.E.S.S. experiments show a much more steeper correlation 
with the index x st 3 (Acharonian et al. 2009). However, in this 
particular case the correlation was observed mostly in a decay 
phase of the flare because of a delay of the Chandra observa- 
tions. Steep correlation means a relatively small change of the 
X-ray flux in comparison with the variation of the gamma-ray 
emission. This means that in an extreme case, where x — > oo 
only gamma-ray activity can be observed. This was reported at 
least two times (e.g. Krawczynski et al. 120041 Blazejowski et 
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al. 120051 1 and is know in the literature as an orphan flare phe- 
nomenon. 

The correlation can be well determined only when an ac- 
tivity event is well observed by two different instruments. This 
requires very good sampling of the recorded light curves (at 
least several observations during a flare) and small error val- 
ues in comparison with the amplitude of the variations. These 
conditions made it possible for only a few cases to precisely 
determine the index of the correlation so far. So was the cor- 
relation sufficiently determined only for three out of nine cases 
analysed by Fossati et al. (I2008I I. Also more recent observations 
of Te V blazars do not give a definitive answer about the correla- 
tion slope (e.g. Albert et al. 2007a Horan et al. 2009 Bonnoli et 
al. 120091 

The observations show that the correlation slope is changing 
form linear to cubic. However, the correlations obtained for a rel- 
atively long period of the observations (weeks or moths) are usu- 
ally linear or slightly more than linear (e.g. Gliozzi et al. 2006 
Albert et al. 2007b), whereas the observations of short flaring 
events (a few hours) give quadratic or even cubic relations (e.g. 
Fossati et al. 120081 Aharonian et al. 2009). Moreover, a scatter of 
the correlated data points seems to be much higher for the long 
period correlations in comparison with the results obtained for 
the short events. The long time correlation contains observations 
of many different flares produced by different components of a 
jet, and this is probably the reason of the scatter. Therefore, the 
short time correlations should provide much better constraints 
for the emission models. Especially interesting are cases where 
the short time correlation is quadratic or more than quadratic. 
Standard one-zone models frequently used to explain the high 
energy emission of blazars are not able to explain such a slope 
of the correlation (Katarzynski et al. [2005 ). In the present work 
we propose a more complex approach, where the emission of at 
least two independent sources is observed at the same time. 

2. The problem of the quadratic correlation 

The most simple model that is able to explain the high energy 
emission of TeV blazars assumes a compact source located in- 
side the jet at a distance of less that 1 pc from the centre. The 
source is filed uniformly by relativistic electrons and a tangled 
magnetic field. The particles spinning around the magnetic field 
lines are producing a synchrotron emission which is usually ob- 
served in the X-ray range. Some fraction of this emission is up- 
scattered to higher energies by the same population of the elec- 
trons. This is the well-known synchrotron self-Compton (here- 
after SSC) radiation that appears in the gamma-ray range. This 
simple scenario was used many times to explain the high energy 
spectra of TeV blazars (e.g. Bloom & Marscher 1996, Ghiselhni 
et al. [19961 Inoue & Tak ahara [19961 Mastichiadis & Kirk [1997l 
Krawczynski et al. 120001 Katarzynski et al. l2001b . 

The intensity of the synchrotron emission is proportional 
to the particle density, whereas the intensity of the SSC radi- 
ation is proportional to the square of the particle density. This 
well-known relationship could explain in principle the observed 
quadratic correlations if only the change of the particle density 
were responsible for the observed activity. The question is how 
realistic such a scenario is. 

An injection of the relativistic particles into the source, 
which increases the density, could in principle explain the 
quadratic correlation during the rising phase of a flare. But the 
particles should be injected simultaneously into the entire vol- 
ume of the source. Moreover, the source volume and the mag- 
netic field inside the source should remain constant during the 



injection phase. The radiative cooling of the particles should also 
be negligible during the injection. All these requirements render 
this scenario unrealistic. We will demonstrate that a more realis- 
tic scenario, which assumes a local injection into an expanding 
source where radiative cooling is important, leads to linear cor- 
relation during the rising phase. 

By analogy systematic energy-independent escape of the 
particles which decrease the density could in principle explain 
the quadratic correlation during the decay phase of a flare. This 
process requires a significantly weaker magnetic field outside 
the source in order to reduce the efficiency of the synchrotron 
emission. On the other hand the radiation field energy density 
outside a spherical (R - radius) homogeneous source at a dis- 
tance of 1/2R above the source surface is only half as weak 
as on the surface (Gould |T9791 >. This means that particles out- 
side the source can still efficiency produce gamma rays through 
the inverse-Compton scattering. In other words, the gamma-ray 
emission will not decay fast enough to produce the quadratic 
correlation during the decay phase. 

The detailed analysis of the correlation for many differ- 
ent scenarios of a single source evolution was performed by 
Katarzynski et al. 2005, where simple analytic formulae were 
derived to describe basic cases, and more complex scenarios 
were analysed through numerical simulations. This analysis 
shows that in all realistic cases we should expect rather linear 
than quadratic correlation. 

Finally, in the case of PK2155-304 the correlation with the 
index x 3 was observed (Aharonian et al. 2009) and this cer- 
tainly cannot be explained by the changes of the particle density 
alone. 

3. Single source vs two sources 

A steep slope of the correlation (x > 2) can be easily obtained 
if we consider simultaneous emission of at least two sources. 
But let us first assume the emission of a single source, where 
the X-ray flux is increasing during a rising phase of a flare as a 
power-law function, where 



= a r 



-at' 



(2) 



is the flux at the beginning (f m ; n ) and the maximum (f max ) of the 
activity respectively. Comparing the two fluxes we can define the 
ratio of the X-ray fluxes in the rising phase 



X,max 'max 



r 



and derive the power-law index 
s — . 

ln(?max) - ln(?min) 



(3) 



(4) 



We can repeat this assumption and calculation for the gamma- 
ray emission, which gives the index for the flux evolution in a 
similar form 



ln(F Te v,i-) 



ln(fmax) - ln(fmin) 



(5) 



Therefore, the slope of the correlation can be described by a sim- 
ple formula 



_ c _ ln(F TeV , r ) 
X ~ s ~ ln(F x , r ) 

that is valid for a single-source emission. 



(6) 
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The formula derived above shows that the correlation slope 
can be increased by an increase of the gamma-ray flux ratio 
(Fxev.r) or by a decrease of the X-ray flux radio (Fx.r)- The first 
option is limited by the nature of the SSC emission which gives 
at maximum a quadratic relation in the specific conditions that 
are rather not realistic, as we discussed already. We decided ac- 
cordingly to test the second possibility and to decrease the X-ray 
flux ratio. This ratio will decrease significantly when we assume 
an additional source of the X-ray emission. In the simplest case 
we can assume a constant level of this additional emission in a 
time F x , const- This gives 

w _ ^X,max ^X, const ^X.max r 

X.' ~ F T~F F = Xr ( ' 

' X,min ' X. const ' X,min 

for Fx, const ^ and shows that Fx q , const must be compara- 
ble or greater than Fx.min to significantly modify the ratio. On 
the other hand, possible gamma-ray emission of this additional 
source should be negligibly small in comparison with the first 
source emission to avoid simultaneous modification of fxev.r- 
Moreover, with additional X-ray emission the correlation is no 
longer a simple power-law function and Eq.|6]can be used only 
for the estimation of the correlation slope. But the derived for- 
mulae explain why the correlation slope can be steeper also in 
more complex situations which are difficult to describe using 
simple analytic relations. 

4. The observations 

In this work we focus on the activity of Mrk 421 observed dur- 
ing the campaign of observations conducted in March 2001 by 
RXTE, HEGRA and the Whipple experiments (Fossati et al. 
2008). The seven days long light curves obtained during this 
campaign were divided into nine periods, and the activity events 
were analysed separately. The best results were obtained for the 
first, fourth and fifth night of the observations (March 18/19, 
21/22, 22/23). 

More than quadratic correlation was observed only during 
the first night. The comparison of the X-ray flux (0.2-10 keV) 
and all the available TeV observations (E > 0.4 TeV) give in 
this case x = 2.26 + 0.25. But the comparison between the X- 
ray data and the TeV observations made only by Whipple dur- 
ing this night gives the almost cubic value of the correlation 
x = 2.84 ± 0.41. What is important, only Whipple was observ- 
ing a strong flare during this night, the light curve obtained by 
HEGRA before the flare shows no activity at all. 

The correlation obtained for the fourth and fifth night was 
less than quadratic but significantly more than linear with x = 
1.56+0.25 and. x = 1.67+0. 16 respectively. Note that the relative 
amplitudes of the flux changes are significantly lower in these 
observations. 

The flare observed during the first night seems to be one of 
the best activity events ever observed in TeV blazars. The evolu- 
tion of the X-ray and gamma-ray flux was well detected during 
the rise and the decay phase of the flare. This gives very impor- 
tant information about the correlation which seems to be very 
similar for the rising and decay phase. Moreover, the excellent 
observations give also information about the spectral evolution 
of the emission. Finally, this activity seems to be quite simple, 
created by one or two sources, whereas the light curves obtained 
during the other nights looks like a superposition of many events 
in time. We consequently selected this particular flare to test our 
models. If a model is able to explain the correlation obtained for 
this extreme flare, it should be able to explain any other correla- 
tion as well. 



5. Injection scenario 

To explain the correlation it is necessary to apply time-dependent 
modeling. Many different models have been proposed to explain 
the evolution of the high energy activity in TeV blazars so far 
(e.g. De rmer[T997l Kirk et al. flt?98l Coppi, & Aharonian [19991 
Kataoka 2000). Most of them assume changes in the particle en- 
ergy spectrum inside a source as the main reason of the activity. 
Some more complex models additionally are taking into account 
the light-crossing time effect (e.g. Chiaberge & Ghisellini 1999 
Sokolov et al. 125041 Graff et al. l2008l >. 

In the first part of this work we use the relatively simple 
model proposed by Chiaberge & Ghisellini ( 1999) and improved 
by Katarzynski et al. (2008). The model assumes particle accel- 
eration by a shock wave inside a jet. The shock is created by two 
or more colliding components of the jet, where some fraction of 
the kinetic energy of the components is converted into random 
energy of the particles. This is the well-known internal shock 
scenario proposed for the first time by Rees (1978 ). However, 
there is no precise description of the acceleration process in the 
model we use. The particle energy increases relatively fast at 
the shock region and after that most of the particles escape into 
the downstream region of the shock. Therefore, the acceleration 
can be approximated as an injection of the particles into some 
volume. In this particular case we assume a power-law energy 
distribution of the injected particles Q(y) = Q m j7 ~", where the 
particle energy is given by E = ym e c 2 . The geometry of the 
shock front is square with the length R and the thickness <K R. 
The duration of the injection process is t — R/c. The injection 
forms a cube like source. This particular geometry was chosen to 
describe the light-crossing time effect. The source volume was 
divided into lOxlOxlOcells. Summing the emissions of the cells 
in the proper way (detailed description in Chiaberge & Ghisellini 
1999) makes it possible to simulate this effect. 

As we already discussed, a single zone model is not able 
to explain a quadratic or a more than quadratic correlation. 
Therefore we propose more complex scenario, where the X- 
ray emission is produced mostly by a relatively large source 
(R ~ 10 16 cm). This source has a relatively small particle density 
and is therefore not able to produce efficient gamma-ray emis- 
sion. To explain TeV gamma rays we use another source that 
is compact (R ~ 10 15 cm) and dense, and is able to dominate 
the emission in this energy range. The two sources are com- 
pletely independent and located at different positions inside the 
jet. But by chance the emission of both sources is observed at 
the same time. For simplicity we use the same injection scenario, 
described above, to calculate the evolution of both sources. 

Note that the scenario which assumes the simultaneous emis- 
sion of two or more sources was already successfully used to 
explain the rapid TeV variability observed in PKS 2155-304 
(Aharonian et al. 120071 ). Compact sources (of a size of about 
10 14 cm) which produce TeV emission several times stronger 
than X-ray radiation can explain activity, where the variabil- 
ity time scale is of about a few minutes. The main problem of 
compact and dense sources - absorption of TeV emission due 
to electron-positron par production, is negligible in such an ap- 
proach (Katarzynski et al. |2008l ). 

The results of our modelling are presented in Fig. [1] where 
we show spectra selected for two different time moments, the 
light curves and the obtained correlation. The source size de- 
termines directly the variability time scale. Accordingly, a more 
extended source produces relatively slow variability, whereas the 
compact source explains the flare. The correlation obtained from 
the modelling does not perfectly reproduce the observed rela- 
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Fig. 1. The activity of Mrk 421 observed by Whi pple an d RXTE 
experiments on 19 of March 2001 (Fossati et al.2008) and the 
results of our modeling which assumes the particle injection. 
The upper panel shows spectra obtained for two different time 
moments: before the main flare and at the top of the flare. The 
middle panels show the light curves where the total emission is 
indicated by the thick lines and, the thin lines show only the ex- 
tended source emission. The dots in the middle panels show the 
time moments, where the spectra were obtained. Lower panel 
show the observed correlation and the results of the modeling, 
where the open symbols indicate the decay phase and the thin 
lines show the correlations calculated separately for each source. 



tion. The correlation index is different for the rising and decay 
phase of the activity. This is directly related to the fact that the 
increment of the emission is produced by the injection, whereas 
the decay is caused by the radiative cooling. The two completely 
different physical processes control the evolution during the rise 
and decay of the activity and the correlation has also different 
slopes during these phases. In Fig.Q]we show the total correla- 
tion obtained for both sources simultaneously and the correla- 
tions calculated separately for each source. The correlation cal- 



culated for the single source has a linear slope (x-f sl 1) during 
the rising phase and and an almost square root slope (x^ ^ 0.5) 
in the decay phase. The total correlation shows that the slopes 
are changing significantly when the emission of two sources is 
observed simultaneously (x-j- =; 1 — > 3 and X[ 0.5 — > 2). This 
helps to explain the observed correlation, but the difference be- 
tween Xf and x^ produced within the single source appears also 
in the total correlation. 

The values of the physical parameters used in the modelling 
which are identical for both sources are: Doppler factor 5 = 20, 
Tmin = 1, Tmax = 10 6 , n = 2. The difference appears in the 
source size R — 1.1 x 10 16 and 1.77 X 10 15 cm, the magnetic field 
strength B = 0.05 and 0.08 G and the density of the injected 
particles f m j2inj = 1.7 x 10 4 and 7 x 10 s cirr 3 for the extended 
and the compact source respectively. The number of free param- 
eters shows that the model we use is quite simple, and this was 
the main reason to chose this particular scenario. A detailed de- 
scription of the model is given in Chiaberge & Ghisellini ( 1999) 
and Katarzynski et al. d2008t . 

Note that we have tested also a scenario which precisely de- 
scribes the acceleration process (Katarzynski et al. 2006). This 
particular model gives a different value of the correlation for a 
single source emission. However, in this case x-f is also signifi- 
cantly different from x^ because two different processes (accel- 
eration and cooling) control the rise and decay of a flare. This ap- 
pears to be the main problem for the models which are trying to 
explain the activity by the evolution of the particle energy spec- 
trum. It seems that within all scenarios which assume the particle 
energy evolution, only adiabatic compression and then adiabatic 
expansion of the source could give x-f = x±. But, this approach 
requires a negligible radiative cooling which makes this scenario 
not realistic. 



6. Doppler boosting effect 

Simultaneous emission of two or more sources can explain any 
slope of the correlation where the slope depends on the relative 
value of the X-ray and the gamma-ray emission of the sources. 
However, the observed correlation seems to be very similar for 
the rising and decaying phase of a flare. This appears to be prob- 
lematic for the models that assume different physical processes 
to explain the two phases of the activity. We decided to test a 
very simple scenario that assumes changes of the Doppler factor 
as a main reason of the flux variations. This idea was already 
proposed several times - for example to explain the activity of 
Mrk 501 (Villata & Raiteri [19991 . 

An intrinsically isotropic emission of a source that travels 
with the relativistic velocity V = Be is confined in a beam 
along the velocity vector. The half-opening angle of the beam is 
<p ^ 1 /r in radians, where F = 1 / yl - B 2 is the Lorentz factor. 
This is the well-know beaming effect that amplifies the observed 
emission. For a given angle 9 between the velocity and the direc- 
tion to the observer the amplification of the emission F v = S^F'^ 
is described by the Doppler factor 6 = 1/(F(1 - Bcosd)). 
Therefore a small increase or decrease of 9 may cause signifi- 
cant variations of the observed flux. This may happen when the 
source travels on a helical or quasi-helical trajectory. Such tra- 
jectories of the jets were observed many times in AGNs and may 
cause quasi-periodic activity observed in blazars (e.g. Rani et al. 
2009 Qian et al. 120091 ). Still for the sake of simplicity we do not 
describe the trajectory of the source in our simulations. This will 
require many free parameters and is not necessary to test our as- 
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Fig. 2. The same activity of Mrk421 as in Fig.Q]and the results 
of our modelling which assume the evolution of the Doppler fac- 
tor. Note that in this particular case the correlation for a single 
source is perfectly linear. This is well visible in the bottom part 
of the correlation panel, where only the background source is 
contributing to the total emission. 



sumptions. We simply describe the change of the viewing angle 
by the cosine function 



0(y) = e u 



AO 



cos(j) + 1 



(8) 



where t9 max = 180/(7rr) and AO is the amplitude of the change. 
Parameter y changes linearly from -n to n during a flare, where 
t a describes the duration of the activity in the observer's frame. 

Note that the amplification of the emission is frequency- 
independent and will always give a linear correlation = 
Xi = 1). This is the main advantage of this approach. Moreover, 
the correlation slope does not depend on the evolution of the 
Doppler factor in time. Therefore, the activity can be simulated 
for example by the source acceleration and deceleration, but this 
approach requires special justification. We assume a constant 



value of the source velocity and small changes of the viewing 
angle that give significant variations of the observed flux. Note 
that this requires a relatively small value of < 10° in general. 

To describe the observed spectra we use a very simple model 
of a spherical and homogeneous source with a uniform density 
and magnetic field. The approach used before which assumes 
simultaneous emission of a big and small source is also used here 
to explain the more than linear correlation. We assume the same 
velocity for both sources that gives Y = 10, and # max = 5.73 deg 
and the same amplitude of the viewing angle change AO = 4.73. 
However, the assumed duration of the activity for the big source 
is t a = 18.5 h, whereas for the small one it is t a = 4.2 h in the 
observer's frame. This indicates that the sources are traveling 
along different trajectories. The difference appears also in the 
radius R = 10 16 and 1.45 x 10 15 cm, the magnetic field strength 
B = 0.1 and 0.05 G and the particle density K = 2.5 x 10 3 
and 2 x 10 5 crrT 3 for the big and small source respectively. The 
particle energy distribution inside the big source is assumed to 
be a broken power-law function N(y) = Ky~" [ for y < ybik and 
N(y) = KyH k "'y-" 2 above the break. We use y min = 1, y brk = 

2 x 10 5 , y max = 10 6 , n\ — 2 and «2 = 4 to describe this spectrum. 
The energy spectrum inside the small source is approximated by 
a power-law function, where we use y m [ n = 1, y max = 2 x 10 6 and 
n — 2. The values used here are very similar to the parameters 
used in the injection scenario. Note that the intrinsic emission of 
the sources is constant in time. This means that the sources are 
in equilibrium, where the radiative cooling is fully compensated 
by the acceleration. The probability of such a situation is rather 
low, but this scenario is physically possible. Dominance of the 
acceleration or the cooling will slightly disturb the perfect linear 
correlation provided by the Doppler boosting effect. However, to 
test an ideal case, where for a single source Xf = x±, we assume 
the equilibrium. 

To compare the boosting scenario with the approach pre- 
sented before we simulate the same activity of Mrk 421. The 
result of the simulation is presented in Fig. [2] This particular 
model provides a better fit for the observed spectra and the light 
curves. However, this is a very simple scenario that has more 
free parameters. A single source emission provides a perfectly 
linear correlation which is well visible in the lower part of the 
correlation panel in Fig. [2] where only the big source emission is 
correlated. Simultaneous emission of both sources gives almost 
an cubic correlation. However, the correlation is not the same for 
the rise and decay of the flare. This is the result of the relative 
shift in time between the light curves produced by each source. 
The flare produced by a single source is symmetric in time in 
this particular scenario. But the maximum of the flare produced 
by the small source is delayed by about half an hour in com- 
parison to the maximum of the big source activity. This causes 
the small difference between jtj and x±. This shows an important 
fact: that the correlation produced by many independent sources 
can be diluted simply by the delays between the flares. This may 
be the reason why the correlation was well determined only in 
a few cases. Finally, the calculated correlation is slightly shifted 
in comparison with the observed correlation. This is the result 
of the delay between the observed X-ray and gamma-ray flare. 
Such a delay cannot be simulated in our simple scenario. 

7. Summary 

Single zone models frequently used to explain emission of 
blazars in the TeV range are not able to explain quadratic or 
more than quadratic correlation between the X-ray and gamma- 
ray emission. So we propose a simple solution to this problem 
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which assumes a simultaneous emission of at least two indepen- 
dent sources. In the first approach we use the model that simu- 
lates activity by the particle injection. This rather classical ap- 
proach was successfully applied in many similar models. But 
our calculations show that the simulated correlation has differ- 
ent slopes for the rising and decay phase of the flare. This seems 
to be a general problem for the models that assume different pro- 
cesses (e.g. acceleration and cooling) to explain rise and decay 
of a flare. This lead us to propose an alternative approach, where 
the activity is produced by the Doppler boosting effect. This sce- 
nario provides a perfectly linear correlation for the single source 
emission in the rising and decaying phase of the activity. But we 
still have to simulate simultaneous emission of the two sources 
to obtain quadratic or cubic correlation. The proposed approach 
has several advantages: 

- it can explain any slope of the correlation, 

- in the extreme case it is possible to explain the orphan flares, 

- the approach does not involve a new model of the emission, 
it uses the standard SSC scenario to explain a single source 
radiation, 

- it may explain why the correlation was well determined only 
in a few cases so far, 

- in was already shown that using this approach it is possible 
to explain also the rapid variability. 

The correlation may give important information about the 
region of the jet where the high energy activity is generated. 
However, many further observations are required to fully under- 
stand the nature of the correlation. There are still many open 
questions. For example, the correlation is observed exception- 
ally, is this normal for TeV blazars or is this rather a problem 
with the observation quality? Is there any preferable slope of the 
correlation? How does the correlation look in different sources? 
Answers for the above questions and for many others will require 
precise, simultaneous observations in the X-ray and gamma-ray 
range. 
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